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ABSTRACT 

We explore the outcome of mass transfer via Roche lobe overflow (RLOF) of Mhs ^ O.SIMq 
pure helium burning stars in close binaries with white dwarfs (WDs). The evolution is driven by 
the loss of angular momentum through gravitational wave radiation (GWR), and both stars are 
modeled using Modules for Experiments in Stellar Astrophysics (MESA). The donors have masses of 
Mne = 0.35,0.4, & O.SIMq and accrete onto WDs of mass A/wd from 0.6Mq to 1.26M0. The initial 
orbital periods (Porb) span 20 to 80 minutes. For all cases, the accretion rate onto the WD is below 
the stable helium burning range, leading to accumulation of helium followed by unstable ignition. The 
mass of the convective core in the donors is small enough so that the WD accretes enough helium-rich 
matter to undergo a thermonuclear runaway in the helium shell before any carbon-oxygen enriched 
matter is transferred. The mass of the accumulated helium shell depends on A/wd and the accretion 
rate. We show that for Afne ^ O.4M0 and A/wd O.SMq, the first flash is likely vigorous enough to 
trigger a detonation in the helium layer. These thermonuclear runaways may be observed as either 
faint and fast .la SNe, or, if the carbon in the core is also detonated. Type la SNe. Those that survive 
the first flash and eject mass will have a temporary increase in orbital separation, but GWR drives 
the donor back into contact, resuming mass transfer and triggering several subsequent weaker flashes. 
Subject headings: stars: binaries: close - stars: novae - stars: cataclysmic variables - stars: white 
dwarfs - supernovae: general 


1. INTRODUCTION 

Core burning helium stars in semi-detached close bi¬ 
naries with WDs are one of the progenitor classes of AM 
Canum Venaticorum (AM CVn) stars, along with donors 
that are helium WDs and evolved main sequence stars. 
They are part of a larger class of binary stars called cat¬ 
aclysmic variables, but are further defined by their short 
orbital periods (5-80 minutes) and a complete absence 
of hydrogen. Their importance stems from common en¬ 
velope stages during their formation, their detectability 
as gravitational wave sources, and their roles as possi¬ 
ble projenitors for h e lium novae, SNe . la, and even SNe 
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is realized/encountered when the initial binary configu¬ 
ration leads to so-called case CB mass transfer: the pri¬ 
mary fills its Roche lobe during its AGB phase (C-type), 
and after a common envelope phase, the second ary fills 


its Roche lobe during its RGB phase (B-type) (Iben & ,,, 

Tutukov (1994)). The two common envelope phases re- (|£UJ^,[£ 


lease a large fraction of the initial angular momentum, 
leaving a helium star in a short orbital period with a 
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WD. The binary is then driven closer together due to 
angular momentum loss via gravitational wave radiation 
until the helium star fills its Roche lobe and begins do¬ 
nating helium to the WD. Longer initial orbital periods 
lead to helium stars that complete their burning before 
making contact. 

When the WD has accumulated enough accreted he¬ 
lium, a thermonuclear runaway begins in the helium 
shell. At this point, the W D will likely eject all th e 


mass above the burning layer (Iben & Tutukov (1991)). 
Depending on the mass of the WD and ot the accreted 
helium shell (and of the fraction of metals in the ac- 
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The goal ot our work is to use a realistic mass trans¬ 
fer model to calculate the He shell masses at ignition 
and assess the likelihood of their detonations. The mass 
loss from these events, which, given survival of the ac- 
cretor, occur multiple times for each system, increases 
the orbital separation (at these low mass ratios), causing 
the stars to temporarily lose contact, until GWR drives 
them together again. Therefore, while the initial com- 
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bined mass of the two compact stars may exceed the 
Chandrasekhar limit, this channel will not result in a 
near Chandrasekhar mass WD. 

In ^we describe the creation and subsequent simula¬ 
tion of the helium star models, compare the secular evo¬ 
lution results to those of previous studies, and explain 
the donor stars’ response to mass loss during different 
stages. Then in ^^we include the behavior of the ac¬ 
creting WD and the conditions leading to the first ther¬ 
monuclear runaway in order to determine the strength 
of the flash and the likelihood of a detonation occurring. 
As a c ase study, we simu lated the binary system CD-30° 
11223 (Geier et al. (2013)) up to the first flash. Addition¬ 
ally, we follow through with the simulations assuming the 
first flash does not unbind the accreting WD, and calcu¬ 
late the effect of mass loss from ejection events on the 
binary parameters in ^ and follow up with conclusions 
and discussion in 

2. DONOR EVOLUTION AND MASS TRANSFER RATES 

These types of binary systems are formed when two 
comparable mass stars, with mass ratio less than 2, are 
born in a close binary. The more massive (primary) star 
will leave the MS and evolve through the RGB and AGB 
phases, donating some of its mass to its less massive com¬ 
panion. The primary star then becomes a WD, while 
the secondary, now the more massive of the two, accel¬ 
erates through its MS evolution. Once the secondary 
leaves the MS and expands as it rises up the RGB, it 
overflows its RL and begins unstable mass transfer. The 
resulting common envelope removes the secondary’s hy¬ 
drogen envelope, revealing a low mass helium star. This 
happens soon enough that the WD’s core is still hot 
(log(Tc/K) > 7.2) at an age of «200 Myr. The two 
stars are then close enough for GWR to drive them into 
contact before exhaustion of helium in the secondary’s 
core as lo ng as Pnrh 5, 2 hours after t he final common 
envelope (iTu t ukov fc Fedoroval (|1989|) , i Tutukov fc Yun- 

elson (|l99U|), llben &: 'Tutukov 

TOMl) , [Y7Egels'3(^ ( [2008t ) . - 


Iben & i'utukov 


Uur 0.35 Mq helium star model is created by evolving 
a 3.0 Mq ZAMS star with solar metallicity, using a mix¬ 
ing length of twice the pressure scale height {umi = 2H) 
and the Schwarzschild criterion with no convective over¬ 
shooting, through the main sequence and into hydrogen 
shell burning until the mass fraction of hydrogen at mass 
coordinate Mr = 0.35 Mq has dropped below 10“®. We 
then artificially remove mass from the surface until the 
mass is 0.35 Mq. The star is left to adjust to the mass 
loss until core helium burning luminosity reaches 1 Lq. 
The 0.4 Mq helium star is created the same way starting 
with a 3.5 Mq ZAMS star. We also created a 0.51 Mq 
sdB star by evolving a 1 Mq ZAMS star up to the helium 
core flash, at which point we turn on a very strong RGB 
wind until the mass drops to 0.51 Mq. The core flash 
then proceeds to stable core He burning. 

We set “tau_f actor=100”, which puts the outer cell of 
the model at an optical depth of r = (2/3) x tau_f actor. 
This setting omits the Hell ionization zone and leads to 
a thin surface convection zone and gives a Tgff consis¬ 
tent with sdB star observations. Setting a tau_factor=l 
leads to a deep surface convection zone, and initial mass 
transfer rates experience a runaway and encounter nu¬ 
merical problems. We leave this issue for future studies. 


2.1. Mass Transfer Rates 



Figure 1. Comparison of M calculated b y MESA’s root-finding 
algorithm (solid line) and using equations Q and Q given by 
component masses and /?He given in MESA output (dashed line). 
Arrows show the direction of evolution. From the case with Mjje = 
O.4M0, = O.6M0 at Porb,0 — 40 minutes. 


The binary interac t ions are first comput ed, using MESA 
(Paxton et al. (2011), Paxton et al. (2013)) version 5118, 
by adding a point-rnass companion, i.e. the evolution of 
the companion star is not computed during this calcula¬ 
tion. We assume no rotation or magnetic braking, and 
conservative mass transfer. The only factors that affect 
the binary separation, a, are mass transfer and loss of 
angular momentum via GWR. To compute mass trans¬ 
fer rates, we used the “Ritter” implicit scheme of MESA 
(Paxton et al. 2015 in preparati on), which im plicitly 
computes the prescription given by Ritter] (1988). 


The resulting mass transfer rate, Mne, matches what 
we would expect from the relation of mass and angular 


momentum loss 

rates, 
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where Mne and M^^d are the masses of the helium 
star and WD, re s pectively, and n = d lnRrtp/dlnMnp, 
(Paczyhski (1967), Nelemans et al. (2 001|)). Angular mo¬ 
mentum loss via GWR is given by (Landau & Lifshitz] 
(119^ as 
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Table 1 

Comparison to previous work by Yungelson 


AtHe/AA© 

Mwd/A t© 

^’orb.o (min) 

tc (Myr) 
Yungelson 

MESA 

Pc (min) 
Yungelson 

MESA 

Yungelson 

MESA 

0.35 

0.5 

20 

1.29 

1.4 

15.96 

15.7 

0.977 

0.974 

0.35 

0.5 

40 

15.99 

17.4 

16.24 

16.1 

0.936 

0.923 

0.35 

0.5 

60 

50.80 

51.1 

17.02 

16.9 

0.871 

0.833 

0.4 

0.8 

40 

9.14 

9.2 

20.61 

21.2 

0.933 

0.925 

0.4 

0.8 

60 

30.21 

30.4 

21.67 

22.0 

0.854 

0.824 

0.4 

0.8 

80 

66.95 

67.4 

22.85 

24.1 

0.751 

0.637 


Note. — Mhc and Mwd are the initial masses of the donor and the accretor, Porb,0 is the initial orbital period, tc is the time before 
the donor fills its Roche lobe, Pc is the orbital period at time of contact, and Yc is the central Helium mass fraction at the time of contact. 



Figure 2. The two models shown here are the 0.51 Mq and the 
0.4 Mq donors. The two timescales for each model are the thermal 
and the mass loss timescales. At the beginning of the mass trans¬ 
fer phase, the thermal timescale is much shorter than the mass 
loss timescale, so the donor star adjusts its thermal structure in 
response to mass loss. As mass transfer rates rise and core lu¬ 
minosity falls, the timescales become comparable, and the donor 
responds adiabatically to mass loss, causing thermal timescales to 
rise and mass loss timescales to fall even faster. 


Figure shows the M computed by MESA and the M 
calculated by equations Q and ( 1 ^ against Porb given 
component masses and Pue (for me derivative n) re¬ 
ported in MESA output for a case with Mue = O. 4 M 0 , 
My/u = 0.6Mq at Porb.o = 40 minutes. The derivative 
n = c?lni?He/dlnAfHe was smoothed in certain regions us¬ 
ing a raised cosine window because of numerical noise. 
The donor hrst starts RLOF at Porb ~ 21 minutes and 
A4 quickly rises to « 3 x lO“®M 0 /yr, follows the arrows 
to a period minimum at 10 minutes, and then expands 
to larger Porb with a declining M. The glitch at wIS.S 
minutes results from the exposure of the once convective 
helium burning core. 

In addition to checking M histories from MESA against 
equations M. and ([ ^, we also compare our results to the 
work of Yungelson (|2008 ). Shown in Table fl] are the age, 
Porb, and the core mass fraction of helium at the time of 
contact for six cases. We agree with this prior work. For 


a hnal comparison to Yungelson’s results, we plotted the 
evolution of the thermal (tkh = GM^ /RL) and angular 
momentum loss (tgw = {J/J)\gwr) timescales for the 
case with initial parameters Mho = 0.35 Mq, = 0.5 
Mq, and Porb.o = 20 min and found go od agreement with 
the same plot shown in Figure 3(a) in Yungelson (2008). 


2 . 2 . Evolution of the Helium Star Donors 

Figure shows the evolution of the thermal and mass- 
loss = M/M, related to tgw by equation 

timescales for two cases with initial parameters (Mne, 
Mwd, Porb,o) = (0.51 Mq, 0.76 Mq, 70.5 min), (0.4 
Mq, 1.0 Mq, 40 min). Models with lower donor mass 
climb to higher M faster because they make contact 
at smaller orbital separations, which means that radii 
and luminosities decrease faster and tkh ~ happens 
sooner in the evolution. This timescale crossing happens 
almost immediately after mass transfer turn-on for the 
two smaller mass donors, while the model with initial 
mass Mue = 0.51 Mq has « 0.2 Mq removed before 
Tkh « 

As the nuclear burning luminosity starts to decrease 
due to mass loss, the core temperature slowly decreases, 
allowing core contraction to increase the central density, 
shown in Figure ^ The bump in Pnuc/T in the middle 
panel of Figure|^is caused by the envelope’s abs orption of 
the nuclear burning luminosity due to mass-loss (Savonije 
|et al.| (|1986|)|Yungdson| (|2008|)), explained more in ^ 
Alter tne buinp in the top and middle panel of Figure 5 
n begins to fall, which, as can be seen from equation ([^, 
increases M. At this point the star is being driven out of 
thermal equilibrium as tkh ^ so further mass loss 
causes the central density to dec rease, hence the ben d 
in the Pc — 7c curves in Figure Isl dSavonije et al.l (|1986|), 
Deloye et al. (2007)). 

The increase in M continues until the period mini¬ 
mum is reached, where n = dlnPHe/dlnAfue ~ 0.14 for 
each model. The decreasing central temperature and 
density keep the star evolving along lines of nearly con- 
stant entropy where th e cores are only mildly degenerate 
(jSavonije et al.| (|1986|), |Deloye et al. (|2007)). The core 
electron degeneracy parameter, rjc = fJ-e,c/kBTc, that 
each of the three models considered here reach when their 
masses drop to 0.2 Mq are (Mhs, Pc) = (0.51 Mq, 3.8), 
(0.4 Mq, 3.1), (0.35 Mq, 2.6). The fact that the donors 
are only mildly degenerate means that they will have 
higher mass transfer rates for a given Porb than their 
fully degenerate cou nterparts during the AM CVn stage 
(Deloye et al. (2007)). 
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Figure 3. Evolutionary tracks of the central conditions in the 
donor stars. The Mjje = 0-51 case has a Mwd = 0.76 Mq 
companion with Porb.o = 70.5 min and donor core electron degen¬ 
eracy parameter rjc = 3.8 (see end of sect. |2.2[ | when Mhb = 0.2 
Mq. The Mjje = 0.4 Mq case has a A^wd = ^-0 Mq companion 
with Porb.o = 40 min and donor core electron degeneracy parame¬ 
ter rjc = 3.1 when Mjje = 0.2 Mq. The Mhb = 0.35 Mq case has a 
Mwu = 1.0 Mq companion with Porb 0 = 20 min and donor core 
electron degeneracy parameter r]c = 2.6 when Mjje = 0.2 Mq. 



0.0 0.1 0.2 0.3 0.4 0.5 

Mass (Mq) 


Figure 4. The models here are the same as plotted in Figures]^ 
andl^ The top section shows the derivative n = dlnRug/dlnMug, 
which was smoothed. The red circles in the middle panel mark 
when the core convection zone disappears. The mass at which the 
period minimum is reached for each model are (MiuUiah ^Porb min) 
= (0.51 Mq, 0.21 Mq), (0.4 Mq, 0.24 Mq), and (0.3 M0,’o.28 
Mq), which corresponds to n Ri 0.14 for each case. 


3. ACCRETION ONTO WHITE DWARFS AND THE FIRST 

FLASH 

As mentioned above, enough helium is accreted onto 
the WD for a thermonuclear runaway to occur. A neg¬ 
ligible amount of helium is burned before the runaway 
occurs, as the stable helium burning range for W Ds of 
this mass is lO~^ M| 0 /yr< M < 3 x lO“®M 0 /yr (Pier- 
santi et al. (2014), Brooks et al. 2015 in preparation), 
a full order of magnitude larger than the accretion rates 
achieved by the cases in this study. The outcome may 
either be dynamical or hydrostatic, with the WD slowly 
expanding to fill its Roche lobe. Under the assumption 
that ejecta take the specific angular momentum of the 
accretor, mass loss causes the binary separation to in¬ 
crease, which temporarily shuts off mass transfer until 
GWR brings the component stars back into contact. We 
chose not to explore the possibility that dynamical fric¬ 
tion within expanding e nvelope removes extra angular 
momentum (Shen (20151). Below we show the results of 
“true” binary runs where both stars are evolved simul¬ 
taneously, including the binary parameters, using MESA 
version 6596. The WD models were created 



log[Orbital Period] (minutes) 

Figure 5. Black solid lines are systems with Mwd = O.SMq, 
Grey solid lines are systems with = I.OMq. The black 

dashed line is from a system with M^d = 1.26M0. The black 
dotted line is from a system with Mwd = 0.6Mq. The grey dotted 
line is from a system with Mwd = 0.5 M(d. The blue line is from 
the system modelling CD-30° 11223, see §3.2| The symbol marking 
the first flash correspond with those in Figure 


in MES A by the same method described in |Wolf et al.| 
(|2013|). We age the WD models to 200 Myr, the relevant 
age tor the 0.35 and 0.4 Mq helium star donors. These 
WDs have core temperatures of a fewxlO^K at the start 
of mass transfer. The WDs with sdB star companions 
could be even colder. 

In this section we show results from cases with (Mne, 
Mwd) = (O. 35 M 0 , O. 8 M 0 ), (O. 35 M 0 , I.OM 0 ), (O. 4 M 0 , 
O. 8 M 0 ), (O. 4 M 0 , I.OM 0 ), (O. 4 M 0 , I. 26 M 0 ), (O.46M0, 
O. 5 M 0 ), (O. 46 M 0 , O. 6 M 0 ) all with initial orbital peri- 
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ods of Porb.o = 40 minutes, and a case with (O.SIMq, 
0.76Mq) and Porb,o = 70.5 to model the observed sys¬ 


tem CD-30° 11223. 
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Figure 6. Accumulated helium shell mass for the first flash. Line 
Types, colors, and markers are the same as for Figure]^ 
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Figure 7. First burst from all “true” binary runs. Blue markers 
are from I.OMq WDs, red markers from O.SMq WDs, green from 
1.26M0 WDs, purple from O.SMq WDs, and cyan from O.GMq 
WDs. Circles are from O.SSMq donors, triangles from O.4M0 
donors, squares from O.46M0 sdB star donors. Filled markers had 
CpT/tnuc < TTf'/C'sound ■ The Orange star-shape is from the case 
study CD-30° 11223. The last (larger) burst on the 1.26M0 WD 
is included as a special case and is discussed in ij^ 


Our immediate goal is to calculate the thermal prop¬ 
erties of the accumulating helium shell on the WD and 


find the shell mass at the time of the thermonuclear run- 


away . Mos t earlier work (see Figure 1 of Shen & Bild- 
sten (2009)) assumed a constant accretion rate, whereas 


these binaries undergo a factor of about three change in 
M during the accumulation phase. Figure shows the 
changing M as the orbital period shrinks for our distinct 
cases. We terminate this plot at the moment of the first 
helium shell thermonuclear runaway. The resulting he¬ 
lium shell masses at the time of ignition (shown in Figure 
\m are not that different than earlier works, but are now 
aerived self-consistently. At these accretion rates, there 
is adequate time for thermal contact with the colder core, 
resulting in a temperature inversion in the helium shell 
and an ignition location above the base of the freshly ac¬ 
creted material. Hence, the convective shell masses will 
always be less than the accumulated masses and we do 
not concern ourselves with possible mixing durin g the 
burnin g into the underlying CO (or ONeMg) core. Piro| 
(2015) recently discussed mixing during accumulation, 
which also looks unlikely at these lower accretion rates. 
Hence, we assume our shells have the abundance distri¬ 
bution of the material from the surface of the donor. 

To calculate the nuclear burning rates we used MESA’s 
co_burn network, which includes alpha chain reactions 
up to ^®Si. The qualitative nature of all bursts ar e qu ite 
similar, so we study one typical case in detail in §3.1| 

In Figure we plot the core mass, which includes all 
the helium below the burning zone, against the con¬ 
vective env e lope m ass. The black lines, from |Shen| 
& Bildsten| (12009 ), represent where the local heating 
timescale [Cpl'/ enuc) equals the (dashed line) local dy¬ 
namical time (Fl/Csound), and (dotted) 10 x the local 
dynamical time, included for uncertainty in detonation 
requirements. The filled markers met the condition that 
the burning timescale became shorter than the sound 
travel time across a hemisphere (cpT/cnuc < Trr/Csound), 
meaning that the burning was at least non-spherically 
symmetric. 


3.1. Typical First Helium Thermonuclear Instability 

The explosion studied here is the first from the case 
with Mne = O.35M0 , Mwd = I.OMq (filled blue cir¬ 
cle in Figure]^. Because the compressional heat leaks 
both out of the surface and into the core, the maximum 
temperature and burning zones are located within the he¬ 
lium shells, instead of the core-envelope boundary. This 
“off-center” is 


constant M { 

Nomoto 

(1982b|, 

Iben & Tutukov 

(1991), 

Woosfey & Kasen 

(2011^. The 

base of the helium shelf 


(where helium mass fraction Y < 0.05), shown by the 
black markers in Figure!^ therefore, lies below the igni¬ 
tion location, and stays fairly degenerate. 

The evolution of the instability proceeds most rapidly 
near the minimum local heating timescale (rheat = 
CpT/cnuc)- Shown in Figure]^ at the base of the con¬ 
vective envelope (coincident with the burning zone) the 
temperature increases and the pressure decreases by a 
factor of «2, and the carbon mass fraction (A 12 ) rises 
from 0.02 to 0.08. The heating timescale, Theat, reaches 
a minimum of ^3.5 seconds, compared to a dynamical 
time (rdyn = H/Csound) of ^0.2 seconds, and the sound 
travel time around a hemisphere, Trr/Csound ~ 6 seconds. 
Convection extends almost all the way to the surface, and 
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becomes fairly inefficient with convective speeds reach¬ 
ing «15% of the sound speed and eddy-turnover times 
(Tconv = H/vconv) drop below 1 second, faster than the 
heating timescale, Theat- This particular flash does not 
become dynamical, as dynamical flash models could not 
be evolved passed a Theat minimum. 



log P (g/cm^) 

Figure 8. Temperature-density profiles from the specific helium 
explosion studied in this section. The 0.0436 Mq of helium sits on 
a I.OMq C/O core. The size of the convective envelope (all the 
mass above the burning zone) is O.OITIMq. The first profile is 
dashed (1530 yrs before min. riieatli the second profile is dotted 
(0.5 yrs before min. Theat)) a-nd the third, fourth and fifth are 
labelled above by time since min. Theat- The black markers show 
the base of the helium shell in each profile (helium mass fraction 
Y < 0.05). The light blue dotted line marks the degeneracy border 
where fie = kgT. 


10 we plot a composition profile of the helium 
Tie accreted material is almost pure helium 


In Figure ' 
envelope. 

with « 1 % and the pre-dynamical convective burn¬ 
ing produces « 2 % and negligible 
As the results of Figures and show, these helium 
shells evolve rapidly and we certainly should not de¬ 
pend on a hydrostatic code to resolve any transition to a 
more dynamical outcome such as a deflagration or det¬ 
onation. However, we have shown that many of these 
first flashes (especially those with Mhb > 0.4Mq and 
Mwd > O. 8 M 0 ) have an adequately massive convecting 
helium shell to trigger dynamical burning (i.e. deflagra¬ 
tion o r detonation). In light of more recent work, specif¬ 
ically [SherT'^F^Iooi^ (|2014|), who showed that using a 
larger nuclear reaction network that allows the proton- 
catalyzed a-capture ^^C(p, 7 )^^N(a,p )^®0 significantly 
reduces that effective "‘He lifetime, detonations are able 
to begin and propagate in much lower mass shells. This 
and other reaction rates depend on the relative abun¬ 
dances in the burning layer. Therefore we need to know 
not only the composition of the accreted matter, but 
also the build-up of burning products from pre-dynamical 
convective burning. 
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Figure 9. Same model as Figure Shown are the evolution of 
four quantities defined at the maximum burning zone, with t = 0 
corresponding to the minimum Theat- The top row shows the local 
heating timescale, Theat = CpTjenuc- The second row shows the 
local pressure, Pb 21 = P^/lO^^erg- cm“^. The third row shows 
the local temperature, T 5 g = T^/IO^K. The bottom row shows the 
local mass fraction of X 12 



Figure 10. Composition and temperature profile of the helium 
envelope for the first burst of the system with Mue = O. 35 M 0 , 
MwD = 1 .OM 0 near minimum Theat 


With the composition information from Figure |10[ 
which, as we stated above, is typical of the bursts stud¬ 
ied here, we can plot each of these bursts on Figure 10 of 
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Shen & Moore (20141, which tells that all the bursts stud- 
ied here will be able to laterally propagate helium deto¬ 
nations. Furthermore, it tells us that the largest bursts 
from the O.SMq and the 1.0Mq WD accretors, and all 
the burst from the 1.26M0 WD accretor, will produce a 
signihcant fraction of isotopes that are radioactive on rel¬ 
evant timescales (X48Cr + -^52Fe + -^56Ni > 0.2). The rest 
of the smaller bursts will be rich in ^®Si and ^^Ga, but 
will mainly consist of unburnt helium, which will likely 
be neutral and unobservable. 


3.2. Case Study: CD-SCP11223 

The system CD-30° 11223 is the only known sdB-fWD 
binary system expected to make contact within the sdB’s 
core He-burning l ifetime. We used data measured by 
Geier et al. (2013) to model this system, which includes 
a U.51A1 q sdB star and a 0.76 Mq WD in an orbital pe¬ 
riod of 70.5 minutes. Comparison of measurements to 
simulated sdB evolutionary tracks in the Teff — log g dia¬ 
gram suggest that the sdB had just recently been formed 
and started the core He-burning phase. 

Given the longer initial orbital period of this system 
compared to the systems in this study, the sdB had a 
longer time for core He-burning to inject heat into the 
envelope, leaving the star in a higher entropy state when 
mass transfer starts. This leads to a lower mass trans¬ 
fer rate during the plateau phase by a factor of a few. 
The lower mass transfer rate allows a larger helium shell 
to build up on the WD before thermonuclear runaway 
occurs. 
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Figure 11. Mass transfer rate shown in black solid line, left y-axis. 
Mass of the WD accretor in red dashed line and mass coordinate 
of the ignition location as a red marker, right y-axis. 


Our results from modeling this system show (Figure 
11) that the helium star fills its Roche lobe after 37.3 Myr 
at an orbital period of 32 minutes. The donor transfers 
O.I 8 M 0 of helium to the WD before thermonuclear run¬ 
away after 46.1 Myr at an orbital period of 16 . 5 min utes. 
Compare this to the prediction of Geier et al. (2013|) of a 
helium shell of 0.1 M 0 that explodes after only 42 Myr at 
an orbital period of 27.4 minutes. The initial mass trans¬ 


fer rat es we computed agree well with those in |Geier et al.| 
(2013), but our model climbs to higher accretion rates by 
about 50% by the first burst. 

Our model shows a O.O 9 M 0 convective envelope on top 
of a (O.76M0 C-0 + O.O 9 M 0 cold He) O.85M0 core, 
which will detonate and produce a significant mass of 
radioactive isotopes. 

4. LATER FLASHES AND SUBSEQUENT EVOLUTION 

If we assume that the system survives the first flash, 
the loss of mass via the ejection of the helium envelope 
will cause the binary separation to increase. GWR will 
then bring the component stars back into contact, and 
the cycle of helium accumulation, ignition, and ejection 
continues until the mass transfer rates fall to such low 


values that another ignition never occurs (Shen & Bild- 
sten| (2009)). To model this, each time helium burning 
on tne W D begins to run away, the helium envelope is re¬ 
moved from the system, taking with it the specihc orbital 
angular momentum of the WD. 


I I I I I I I I I I I I I I I I I I I I I. 



6 8 10 12 14 16 18 

Myr 


Figure 12. Four true binary runs all started at Porb o = min¬ 
utes. The WD accretors experience between 3 and 9 bursts, each of 
which temporarily shut-off mass transfer. Mass transfer rates are 
shown by the black solid lines and the orbital periods are shown 
by the red dashed lines. All panels have the same scale. 


We show the effects on the mass transfer and orbital 
period histories c aus ed by including the mass ejection 
episodes in FigurefT^ The 0.35Mq donor models lack the 
plateau feature in toe O. 4 M 0 donor models because the 
thermal timescale (tkh) drops below the mass transfer 
timescale ("ry^) much sooner for lower mass donors. The 
higher mass WDs have lower envelope masses (AMenv) 
necessary for helium ignition, so these models have more 
bursts. This behavior can also be seen in Figure |13[ 
which shows the mass of helium shell in the black lines, 
and the luminosities of the helium star and WD in red 
and blue, respectively. 

Before contact, the helium star is more luminous than 



































































the WD. After contact, that relation is reversed, and 
Awd IOLq. This approximately matches the com- 
pression luminosity, T ^nmn ~ SksTciM)/fimp, given in 


Bildsten et al. (2006). This does not include the accre¬ 


tion luminosity, Lacc ~ GMwDAf/i^wD ~ lOOL©, which 
is released in the accretion disk and will outshine both 
stars. 



Figure 13. Same four binary runs as in Figure [T^ Mass of the 
helium envelope on the WD (AMgnv) shown in black lines, lumi¬ 
nosities of the helium star and WD in red and blue, respectively. 


For bursts that remove enough mass and mass transfer 
is shut off for long enough, the WD can be seen to dim to 
the original cooling curve, before climbing back up to the 
compression luminosity once mass transfer begins again. 
During these same periods, after a burst but before re¬ 
contact, the donor luminosity increases. This happens 
because during mass loss periods, the nuclear burning 
lumino sity in the donors c ore is absorbed in the stellar en¬ 
velope (|Yungelson| (|2008| ). In the donor’s envelope there 
is no nuclear burning, so the luminosity’s slope is domi¬ 
nated by mass loss [dL/dm « {M/AM){kBT/mp) < 0, 
where AM is overlying mass). When mass loss sud¬ 
denly shuts off, this mass loss term disappears giving 
dL/dm « 0 in the envelope, meaning the trapped heat 
gets released on the timescale ^ 10 ^ years. 

The features seen in the bottom two panels after 16 
Myr result from the exposure of the once convective he¬ 
lium burning core. When sharp compositional gradients 
are exposed, the mass transfer rates experience a spike. 

For all cases studied here, except the one with Mwd = 
1 . 26 M 0 , the first burst is the largest because the helium 
is accumulated at the lowest rate. As the mass transfer 
rates decrease after the period minimum, they quickly 


become too low to trigger helium ignitions (Shen & Bild- 
(2009)) and the WD ends up with a cold, thick he- 


sten 


lium shell. The case with the Mwd = 1.26M0 accretor is 
the only case in which the last explosion is the strongest 


(Figure]^, as the ignition mass is much lower for the 
more massive WD accretors. In fact, the last explosion 
on the Mwd = 1 . 26 M 0 accretor is the most powerful one 
studied here. After this last explosion, likely as a .la SN, 
the donor does not have enough mass left to increase the 
mass of the accretor to the Chandrasekhar mass, so an 
accretion induced collapse is not expected through this 
channel. 


5. CONCLUSIONS 

We present the first self-consistent calculation of mass 
transfer between a low-mass helium star and a WD in a 
short-period binary and the resulting helium shells and 
their explosions on the surface of the WD. We started by 
simulating the secular evolution of low mass helium stars, 
that is, evolved with a point-mass binary companion, to 
see how structure and mass-loss feed back on each other. 


We found, confirming Yungelson (2008), that during the 
first stage of mass transfer, the donor maintains ther¬ 
mal equilibrium, as tkh < Tm, leading to a relatively 
constant M during the so-called plateau phase. Even¬ 
tually the donor loses enough mass to effectively shut¬ 
off nuclear burning and reverse the timescale inequality, 
meaning the star responds adiabatically to mass-loss. At 
this point mass transfer rates rise until mass transfer ef¬ 
fects become dominant over GWR in determining the 
orbital separation, and the system reaches a period min¬ 
imum of 9-10 minutes, after which the binary separation 
increases. Core evolution proceeds adiabatically along 
nearly constant entropy trajectories that are mildly de¬ 
generate until nearly all the mass of the donor has been 
removed. 

When we include the evolution of the accreting WD 
in the binary simulations, helium shells build up on the 
WD until a thermonuclear runaway develops. In study¬ 
ing these models, we find that the first flash is usu¬ 
ally the largest, and the first flashes from systems with 
Mne > O. 4 M 0 and Mwd > I.OM 0 meet the require¬ 
ments for triggering dynamical burning and developing a 
de flagration or detona tion. Using the conditions derived 
by Moore et al. (2013) and Shen & Moore (2014), all the 
flashes studied here, including the “later' hashes, are 
massive enough to sustain a laterally propagating deto¬ 
nation within the He shell. Furthermore, the first bursts 
from the O.SMq and the I.OM 0 WD accretors, and all 
the burst from the I. 26 M 0 WD accretor, will produce 
a significant fraction of isotopes that are radioactive on 
relevant timescales. 

We also follow through with some simulations assum¬ 
ing that the first helium flash does not unbind the WD. 
The ejected mass increases the binary separation, tem¬ 
porarily shutting off mass transfer until GWR brings the 
stars back into contact. This series of events repeats it¬ 
self until the stars are far enough apart and mass transfer 
rates are low enough to prevent runaways in the helium 
shells. 

We thank Kevin Moore for helpful discusions regarding 
helium detonations, along with Dean Townsley and Ken 
Shen for their useful comments. We would also like to 
thank the referee for their many helpful comments. This 
work was supported by the National Science Foundation 
under grants PHY 11-25915, AST 11-09174, and AST 12- 
05574. Most of the simulations for this work were made 
possible by the Triton Resource, a high-performance re- 
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percomputer Center at UC San Diego. 
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